Ferromagnetic insulating phase in Pri_ x Ca x .Mn03 
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A ferromagnetic insulating (FM-I) state in Pro.7sCao.25Mn03 has been studied by neutron scat- 
tering experiment and theoretical calculation. The insulating behavior is robust against an external 
magnetic field, and is ascribed to neither the phase separation between a ferromagnetic metallic 
(FM-M) phase and a non-ferromagnetic insulating one, nor the charge ordering. We found that 
the Jahn- Teller type lattice distortion is much weaker than PrMnOs and the magnetic interaction 
is almost isotropic. These features resembles the ferromagnetic metallic state of manganites, but 
the spin exchange interaction J is much reduced compared to the FM-M state. The theoretical 
calculation based on the staggered type orbital order well reproduces several features of the spin 
and orbital state in the FM-I phase. 

PACS numbers: 75.47.Lx, 75.47.Gk, 75.25.+Z, 75.30.Ds 



I. INTRODUCTION 

Since the discovery of the high-T c superconductivity in 
layered cupurates, doped Mott insulators with perovskite 
structures become one of the major targets in recent stud- 
ies of condensed matter physics. Hole doped perovskite 
manganites, J?i_ XJ A x Mn03 (i? is a trivalent rare-earth 
ion and A is a divalent alkaline-earth ion) , have been ex- 
tensively studied since the discovery of the colossal mag- 
nctorcsistancc (CMR) effect. Recent studies devoted to 
clarify an origin of the CMR effect have revealed that 
the physics of manganites is controlled by an interplay 
between spin, charge and orbital degrees of freedom. 

At x = in Ri- x A x MnOs, a Mn ion shows an elec- 
tron configuration of t^g^g where the e g orbital degree 
of freedom is active. Actually, the mother compound 
_RMn03 shows an alternate ordering of the d^ x 2_ r 2 and 
d 3y 2_ r 2 orbitals associated with the cooperative Jahn- 
Teller distortion in the ab plane of the orthorhombic 
Pbnm lattice. The layer-type (^4-type) antiferromagnetic 
(AF) spin ordering is realized at low temperature by the 
anisotropic supcrcxchangc interaction under the orbital 
ordering (QO) ^aa This anisotropy was experimentally 
confirmed in LaMnC>3 by neutron scattering studiesA^ 
By substituting R by A in i?MnC>3, holes are introduced 
in the system. In the simple double exchange scenario, a 
kinetic energy gain of the doped holes brings about the 
ferromagnetic metallic (FM-M) ground statei^ It is also 
supposed that the 00 is destroyed and the orbital liq- 
uid/disordered states are realized in the FM-M phased 
The isotropic characters in the crystal latticed and the 



spin exchange interactions^*^ realized in this phase may 
be attributed to this orbital state. 

However, in the actual compounds, the doping depen- 
dence of the electronic ground state is not so naive. Be- 
tween the ^4-type AF phase and the FM-M phase, the 
ground state is ferromagnetic but insulating* The fer- 
romagnetic insulating (FM-I) state can not be explained 
by the simple double exchange picture. Although such a 
contradiction remains to be solved, the FM-I phase has 
not been much studied so far. This paper provides a de- 
tailed information about the character of the FM-I phase 
by a neutron scattering study and a theoretical calcula- 
tion on Pri-^Ca^MnOs, which is a typical example for 
a FM-I phase with a relatively wide hole concentration 
region O.f 5 < x < 0.3iii 

One of the interpretations of the FM-I phase is based 
on a phase separation between the FM-M phase and a 
non-FM insulating phase. In a recent trend of the re- 
search of CMR manganites, the phase separation phe- 
nomenon is regarded as an important factor for the large 
resistivity change accompanied by the insulator-metal (I- 
M) transition. A plausible origin of the non-FM insu- 
lating phase is the CE-type charge ordering (CO). This 
is a checkerboard-type CO of Mn 3+ and Mn 4+ accom- 
panied by a zig-zag ordering of d^ x 2_ r 2 and d 3y 2_ r 2 or- 
bitals on Mn 3+ sites and the CE-type AF spin ordering. 
The C7?-type CO is most stabilized at x = 1/2, but it 
is exceedingly robust against hole doping. Actually, in 
Pri-zCa^MnOs, the Ci?-type CO is well developed for 
x > 0.3 adjacent to the FM-I phase, and coexists with 
the FM component liiiiSii 3 . Recently, Dai et al. showed 
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by neutron scattering measurement that the short-range 
order of the CE-type CO exists in the FM-I phase of 
Lai-zCa^MnOsjA^ though the detail relation between 
these charge correlations and the transport is not ob- 
vious. 

The phase separation picture requires substantial vol- 
ume fraction of non-FM insulating phase to make the 
compound an insulator overcoming the metallic con- 
duction in the FM phase. Uehara et al. stud- 
ied the phase separation between the FM-M phase 
and the charge-ordered insulating (CO-I) phase in 
La5/8_j / Pr y Ca3/8Mn03. They showed that the concen- 
tration of the FM-M phase at the FM transition, which 
was evaluated from the ratio of the saturation mag- 
netic moment to the full moment of a Mn ion, agrees 
with the three-dimensional (3D) percolation threshold 
10-25%^ However, the FM moment of -2 [i B in 
Pro.7Cao,3MnO j 12 i 13 i 16 is too large to account for its in- 
sulating character. On the other hand, there is another 
possibility that the insulating nature does not arise from 
the non-FM phase but it is intrinsic to the FM phase. 
Recently several kinds of COs related to the hole concen- 
tration of a; = 1/4 were theoretically proposed. Hotta 
et al. claimed that Pr 3 /4Cai/4Mn03 becomes a FM-I 
state by a CO consisting of a checkerboard charge or- 
dered plane and a purely Mn 3+ plane stacking alternately 
along the z axistf As a result, the unit cell of the CO 
becomes \[2 x \/2 x 4 in terms of the cubic perovskite 
unit (We will refer this type of CO as (tt, it, 7r/2)-type 
CO.). Mizokawa et al. proposed that in the lightly 
doped manganitcs there exist "orbital polarons." Here 
Mn 4+ ions are surrounded by ferromagnetically coupled 
nearest neighbor Mn 3+ ions where the e g orbitals are 
pointed to the central Mn 4+ siteiA& They claimed that 
in Pr 3 /4Cai/4Mn03, the orbital polarons form a body- 
centered cubic lattice whose unit cell has a size of 2 x 2 x 2 
in unit of the cubic perovskite lattice, and a FM in- 
sulator is realized^ Neither kind of CO, however, has 
not been detected experimentally so far. Instead, En- 
doh et al. claimed by resonant x-ray scattering (RXS) 
and neutron scattering studies that the FM-I phase of 
La 8 gSi'o.i2Mn03 shows an 00 of a mixture of d 3z i_ r 2 
and d x 2_ y 2 orbitals. They also proposed that this 00 is 
an origin of the FM-I statei^2i2I As for the orbital state of 
the FM-I phase in Pri_ x Caa;Mn03, Zimmcrmann et al. 
advocated by a RXS study that the same 00 as that in 
LaMn03 (LaMn0 3 -type 00) develops below about 500 
K in Pro.75Cao.25Mn03i2£ In contrast to the FM spin or- 
dering in Pr .7 5 Cao.25Mn03, however, LaMn03-type 00 
is generally accompanied with the ^4-type AF spin order- 
ing. 

The complicated status of the current research on the 
FM-I state of manganites reminds us the necessity of a 
comprehensive study on spin, charge, and orbital states 
of the FM-I phase of manganites. For this purpose, neu- 
tron scattering is a useful tool to study CO and 00 be- 
cause of its sensitivity to the lattice distortions induced 
by CO and 00. In the present study, we performed clas- 



tic neutron scattering measurements on a single crystal 
of Pro.75Cao.25Mn0 3 to elucidate the superlattice peaks 
and the lattice distortions due to CO and 00. The 00 
affects the exchange interactions between spins due to 
the anisotropic character of the e g orbitals. Therefore, 
through the estimation of the exchange interactions by 
measuring the magnon dispersion relation, we can char- 
acterize the orbital state. So, we also performed inelastic 
neutron scattering measurements to determine the spin 
exchange through the measurements of the spin wave ex- 
citations. The theoretical analyses provide a unified pic- 
ture for the spin and orbital states in the AF, FM-I and 
FM-M phases. 

The rest of this paper is organized as follows. In the 
next section the experimental procedure is described. 
The results and discussions are described in Sec. Ill, 
where the property of the magnetic and crystal structures 
in Pro.75Cao.25Mn03 are discussed. In Sec. Ill- A, we will 
argue the possibility of the phase separation, and in Sec. 
III-B we will characterize the FM-I state by measure- 
ments of the structural properties and spin fluctuations. 
Section III-C is devoted to the theoretical investigation 
for the magnetic interaction and the orbital state in mod- 
erately doped manganitcs. A summary of the work is 
given in Sec. V. 



II. EXPERIMENTAL PROCEDURES 

Single crystals of Pro.75Cao.25Mn0 3 and PrMn03 were 
grown by the floating zone method. A stoichiometric 
mixture of Pr 6 0n, CaC03, and Mn 3 04 was ground and 
calcined at 1300 °C in air for 12 h. Then the resulting 
powder was pressed into rods and fired again at 1300 °C 
in air for 12 h. The crystal growth was performed in 
a floating-zone furnace in a flow of oxygen gas. We 
precharacterized the samples by a powder x-ray diffrac- 
tion measurement and confirmed that the samples are 
single phase. Then the samples were cut into the size of 
6 mm</> x 30 mm for the present neutron scattering study. 
The transport data were taken using Quantum Design 
PPMS on a sample cut from the same batch as used in 
the neutron study. 

The neutron scattering experiments were performed 
using triple axis spectrometers GPTAS and HQR of In- 
stitute for Solid State Physics, University of Tokyo, in- 
stalled at the JRR-3M research reactor in JAERI, Tokai, 
Japan. The most of the measurements were done at GP- 
TAS, while HQR was utilized when a high Q resolution 
was required. At both spectrometers, the 002 reflection 
of pyrolytic graphite (PG) was used for the monochroma- 
tor and analyzer. At GPTAS, we selected a neutron wave 
length of kf = 3.81 A -1 . For elastic scattering measure- 
ments, a combination of horizontal collimators of 20'-40'- 
40'-80' (from reactor to detector) was adopted and two 
PG filters were placed before the monochromator and 
after the sample to suppress contaminations of higher- 
order harmonics. For inelastic scattering measurements, 



3 



o 

CD 
(/) 
O 
LT) 



C= 
=3 
O 

o 



CD 



a 



CD 
O 

-t— ' 

w 

'w 1 
CD 1 
DC 

1 



Pr .75Ca .25 Mn O3 



1000 



Pr a75 Cao.25Mn03 



10° 



1 


| i i i i | i i i 

(a) : 


\ T 


(200) 

C 


" \ > 


t 








i \'\ I I 1 I I 
- '*\ 


i i i i i i i i i i : 


\\\ 

r A 


(b) : 


\ \\ 




- '> v\ 


-0T : 


\ Ov 


- - 2T ! 


\ \ ^ 
! \ \ 


7T : 

- N \ 







50 100 150 200 
Temperature (K) 



250 



FIG. 1: (a) Temperature dependence of the scattering in- 
tensity of the 200 peak, (b) Temperature dependence of the 
resistance under the magnetic field of H = T, 2 T, and 7 T. 



we chose 40'-40'-40'-40' combination of collimators and 
set a PG filter after the sample. At HQR, which is on the 
thermal guide tube, neutrons of wave length hi = 2.56 
A -1 and horizontal collimation of blank-40'-80' (from 
monochromator to sample) were used together with a PG 
filter before the sample. The crystals were mounted in 
aluminum cans filled with helium gas. The temperature 
of the samples were controlled by conventional and high- 
temperature-type closed cycle He-gas refrigerators. The 
measurements were performed in the (h, 0, 1) and (h, h, I) 
zones of reciprocal space of the Pbnm lattice. 



III. RESULTS AND DISCUSSIONS 

A. Possibility of the phase separation 

First, let us characterize the sample by showing the 
results of the temperature (T) dependences of the or- 
der parameter of the magnetic order and the resistance. 
Figure ^a) is the T dependence of the scattering inten- 
sity of the fundamental 200 peak. The intensity starts 
to increase below the Curie temperature Tc ~ 130 K 
due to the FM spin ordering, consistent with a previous 
study^i To determine the critical exponent, we fit the 
data to a power law, (1 — T/Tc) 213 , varying the lower 
limit of the temperature range of the fitting, Xi; m . (3 was 
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FIG. 2: Temperature dependence of the intensity of the CE- 
type charge/orbital order peak at Q = (3.275,3.275,3.55). 
The intensity due to the Huang scattering measured at Q = 
(3.33, 3.33, 3.66) is subtracted. 



determined to be 0.36(1) as 7]; 



T c = 129.5(5) K. 



This value of the critical exponent is very close to that 
for the 3D FM Heisenberg model, 0.367, and similar to 
that for the FM-M manganitcs such as Lao^Sro^MnO^^ 
and Lao.7Sro.3Mn03i2iS The magnetic moment was de- 
termined from a powder sample which was prepared 
by crushing a melt-grown crystal. The FM moment is 
3.32(3) hb per Mn site at 15 K, oriented parallel to the 
b axis. This value of the FM moment is nearly equal to 
the formula value of 3.75 hb- This means that even if a 
non-FM region existed in the sample, its volume would 
be extremely small. 

Figure Q^b) shows the T dependence of the resistance 
at several magnetic fields. For H = T, the sample is in- 
sulating for all temperature range measured, and a small 
cusp is observed at Tc- By applying the magnetic field, 
this cusp vanishes, and the resistance is reduced from 
that at H = T. However the insulating behavior still 
persists under the magnetic field up to H = 7 T. The ro- 
bustness of the insulating state contrasts markedly with 
the CE-type charge ordered phase for x > 0.3, where the 
CO melts by a magnetic field of a few teslasiii In the 
present case, because the spins are already fully aligned 
ferromagnetically at H = T, an external field cannot 
make a significant change of the transports. 

As mentioned above, the magnetic and transport 
properties suggest that the insulating behavior of 
P r o.75Cao.25Mn03 does not originate from the non-FM 
region, but it is rather intrinsic to the FM state. To 
confirm this point, we performed the neutron diffrac- 
tion measurements and surveyed reciprocal points to de- 
tect signals due to the CE-type CO, and found the lat- 
tice superlattice peaks due to the CE- type CO. Figure 
121 shows a T dependence of the scattering intensity at 
Q = (3.25,3.25,3.5), which corresponds to a CE-type 
CO peak (5, 1.5, 0)2^ of a different domain of the sample. 
Because diffuse scattering due to single polarons (Huang 
scattering) overlaps at this position^ the intensity of 
the diffuse scattering estimated at Q = (3.33,3.33,3.66) 
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FIG. 3: Profiles of the 220/004 doublet along the [110] direc- 
tion at 100 K and 560 K. The data are collected at the HQR 
spectrometer. Solid lines are fits to two Gaussians. 



is subtracted from the data, and intensity of 100 counts 
is added to compensate the difference in background in- 
tensity at the two positions. The CE-type CO develops 
below T ~ 320 K. ft is drastically suppressed below Tq, 
although the sample remains insulating [Fig.^b)]. The 
small cusp of the resistance at Tc may be related to the 
melting of the CZ?-type CO, but the overall insulating 
feature of the resistance, however, cannot be explained 
by the CE-typc CO. This is in contrast to the FM-M 
region of manganitcs, where the disappearance of the 
short-range Ci?-type CO accompanies the transition to 
the FM-M state. 14 i 26 i 27 We think that the existence of 
the CE-type CO is due to a small mixture of a higher Ca 
concentration phase. We also found the pseudo Ci?-type 
AF peaks concomitant with the Ci?-type CO in the re- 
gion of x < 1 /2i 12 i 13 i 16 i 28 Their intensity is of the order of 
10 -3 in comparison with that of the FM peak, indicating 
that the volume of the CE-type CO phase is negligible. 



B. Character of the homogeneous FM-I state 

The results described in the previous section indicates 
the phase separation between the FM-M phase and the 
non-FM insulating phase is not important for the FM-I 
phase of Pri-^Ca^MnOs. The FM-I state is an intrin- 
sic state of the low doped region of manganites. Con- 
sequently, to understand the homogeneous FM-I phase, 
we will characterize the charge and orbital states in this 
phase. 

First, we tried to clarify whether the FM-I state is a 
FM charge ordered state as theoretically proposed. For 
the (tt, 7r, 7r/2)-type CO proposed by Hotta et al.fH the 
superlattice peaks are expected at Q = {h, k, I) with h + 
k = odd and I = half integer. To detect this type of CO, 
we surveyed (m, 0, 1/2) and (5, 0, n/2) positions where m 
and n are odd integer, but could not find any peaks. This 
result is negative for the existence of this type of charge 
order in Pro.75Cao.2sMn03. 



(a) Pro.75Ca .25Mn0 3 
(220)* T = 8K 
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FIG. 4: Profiles of the 220/004 doublets along the [110] di- 
rection for (a) Pro.75Cao.25Mn0 3 and (b) for PrMn0 3 at 8 K. 
The horizontal axes indicate the amplitudes of the scattering 
vector Q. Solid lines are fits to two Gaussians. 



The orbital polaron lattice proposed by Mizokawa 
et ali 18 i 19 will produce superlattice reflections at Q = 
(h, k, I) with (h + k = odd and I = even) or (h, k = 
half integer and I = odd). We found peaks at several po- 
sitions satisfying these conditions and forbidden by the 
Pbnm symmetry. Unfortunately, however, these posi- 
tions coincide with those with integer indices considering 
the twinning of the crystal. Because additional small tilts 
of Mn06 octahedra can produce Bragg peaks at such po- 
sitions, we cannot determine whether the orbital polaron 
lattice exists or not from the present study, though no 
clear relation between the intensities of these peaks and 
the resistance was observed. 

Next, we will investigate the orbital state of the FM- 
I state. According to the recent RXS study, resonant 
signals observed in Pr .75Cao.25Mn03 is consistent to the 
same OO as LaMnOs, i.e. alternate ordering of d 3x 2_ r 2 
and G?3j,2_ r 2 orbitals in the ab planed LaMn03-type OO 
is expected to produce significant anisotropies in crystal 
structure as well as spin exchange: The in-plane lattice 
constants a and b become longer than the out-of-planc 
lattice constant c, and the exchange interactions within 
the ab plane will be stronger than that along the c axis. 
Therefore, measuring these quantities should be useful 
for quantitative understanding of the orbital state. 

Figure shows neutron diffraction profiles around the 
fundamental (2,2,0) position along the [110] direction at 
100 K and 560 K. Because of the twinning, the 004 peak is 
also observed at slightly higher Q than the 220 peak. At 
T = 560 K, both peaks locate so close that the observed 
profile forms almost a single peak. On the other hand, 



5 



Pr 0.75 Ca 0.25 Mn °3 




(0,0,11) fi&2) 



FIG. 5: Profiles of constant-energy scan of the magnon spec- 
tra along the [001] direction (left) and the [110] direction 
(right). The background is subtracted from each data. The 
vertical axes are shifted by 40 counts for each energy. Solid 
lines are fits to the Lorentzian spectral function convoluted 
with the instrumental resolution. 



this peak splits into two peaks at T = 100 K due to small 
increase of the a or b axis and large decrease of the c axis. 
This change of the lattice constants is consistent to the 
cooperative Jahn- Teller distortions accompanied by the 
LaMnC>3-type OO. We measured T dependences of sim- 
ilar profiles and found that the structural change gradu- 
ally develops below ^450 K. This temperature coincides 
with the temperature where the orbital signals measured 
by the RXS study starts to increase^ and corresponds to 
the transition temperature from the O phase to the O' 
phase identified in Ref. |2^. These results qualitatively 
suggest the developing of the LaMn03-type OO below 
T ~ 450 K. However, the amplitude of the lattice distor- 
tion is fairly small compared to the mother compound, 
PrMnOa, where a typical LaMn03-type OO is expected. 
This can easily be assured by comparing the profile of 
the 220/004 reflections for both compounds. In Fig. 0] 
we depicted the profiles of the 220 and 004 reflections of 
Pr .75Cao.25Mn0 3 [Fig. Ufa)] and PrMn0 3 [Fig.Bfb)] at 
T = 8 K. We show the amplitude of the scattering vec- 
tor Q as the horizontal axis for easier comparison of two 
data. 

Figure [S] is typical profiles of constant-energy scans 
of spin wave excitations along the out-of-plane direc- 
tion [001] (left panel) and the in-plane direction [110] 
(right panel) of Pr .75Cao.25Mn03. The data were col- 
lected near the FM Bragg point (0,0,2) at T = 8 K. 
The magnetic origin of these excitations were confirmed 
by measuring their Q dependence. The peak positions 
get away from the zone center as the energy increases 
according to the dispersion relation of the spin wave. 
From these measurements, we obtained the dispersion 
relation of the spin wave along the two directions, which 
is shown in Fig. |SJ Note that the momentum trans- 
fer q = (0, 0, 2n) has almost the same magnitude as 



Pr 75 Ca .25MnO 3 



E 1 







- • (0,0, Tl) 


/ - 


m- 

















0.1 0.2 0.3 0.4 0.5 

2 1.8 1.6 1.4 1.2 1 
(0, 0,11) 

FIG. 6: Spin wave dispersion relations in Pro.75Cao.2sMn03 
along the [110] and [001] directions at 8 K. Solid line is a dis- 
persion curve for the FM Heisenberg model with an isotropic 
nearest neighbor coupling. 



q = (n, n, 0) in the Pbnm notation. In contrast to the 
LaMnOsj** two dispersion curves well coincide, indicat- 
ing that the spin exchange is isotropic (J a b ~ Jc)- This is 
consistent to the critical exponent (3 [Sec. Ill- A], which 
is almost the same as the value for a 3D system. The 
observed dispersion curves can be described by the con- 
ventional FM Heisenberg model: H — — YLij Jij &i ' Sj > 
where Jij is the exchange integral between spins at site 
Hi and Rj. In the linear approximation, the disper- 
sion relation is given by Ticu(q) = A + 25 f [J(0) — J(q)], 
where J(q) = J2j Jij exp[iq ■ (Ri — Rj)] and A accounts 
for small anisotropics. Taking only the isotropic nearest 
neighbor coupling J into account is sufficient to repro- 
duce the observed data, and we get the solid line in Fig.El 
with 2 JS = 4.0(2) meV and A = 0.0(3) meV. 

The isotropic character of the spin system shows an 
interesting resemblance among the orbital-ordered FM- 
I phase and the orbital-disordered FM-M phase. How- 
ever, one can distinguish them by the magnitude of the 
spin exchange. Clear enhancement of the spin stiffness 
D(T) on the transition from the FM-I phase to the FM- 
M phase has been observed in a hole concentration de- 
pendence of Lai-^Sr^MnOj^ and Lai^ajCa^MnOg^S by 
neutron scattering studies. The same phenomenon on 
the field-induced I-M transition in Pro.7Cao.3Mn03 has 
also been observed^ In these materials, D(0) in the in- 
sulating phase is «50 meVA 2 and that in the metallic 
phase is «150 meVA 2 . If we estimate D(0) of the present 
P r o.75Cao.25Mn03 compound from J at 8 K by a rela- 
tion D = 2 J Sa^. where a c is the lattice constant of the 
cubic lattice, D(0) becomes ~60 meVA 2 . This value is 
almost same as that of the insulating Lai-^Sr^MnOs, 
Lai-zCaajMnOs, and Pro.7Cao.3Mn03. This fact evi- 
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dences the common ground of the FM-I phase of the 
manganites. 

The Jahn- Teller type lattice distortion means that 
there may exist a LaMn03-like 00 as reported by the 
RXS study. Its smallness and the isotropy in </, however, 
suggest that the electron distribution is more isotropic 
than that in the LaMn03-type 00. Similar situation 
is realized in the FM-I phase of La .g8Sro,i2Mn0 3! 20 i 32 
where the staggered orbital ordering of (|d3 Z 2_ r 2) ± 
\d x 2_ y 2)) I 'y/2 is proposediSSiSi In the next section, we will 
theoretically investigate the orbital state Pri-^Ca^MnOa 
and construct the spin and orbital phase diagram. 



C. Theory 

Spin and orbital states in manganites with moderate 
hole doping are examined theoretically from the view 
point of strong electron correlation. We adopt the t-J 
type Hamiltonian where the double degeneracy of the e g 
orbitals is taken into account. The model Hamiltonian 
consists of the following four termst^ 

7i = Tit + Ttj + Tin + T^-af- (1) 

The first and second terms are the main terms corre- 
sponding to the t and J terms in the t-J model, respec- 
tively, given by 

Wt= E tt^Jn'v+H.c, (2) 

(ij)jy'a- 



and 




di ia is the annihilation operator of the e g electron at site 
i with orbital 7 and spin a. This operator is defined in 
the restricted Hilbert space where the e g electron num- 
bers at each Mn site are limited to be one or less. The 
prefactor tJJ is the transfer integral between orbital 7 
at site i and 7' at j. Tij describes the superexchange in- 
teractions between the nearest neighboring Mn spins and 
orbitals. The spin and orbital degrees of freedom for the 
e g electron are represented by the spin operator 

7(7(7' 

and the pseudo-spin operator 

Ti = — ^2 ^I 7 o-^77' C '»7'ct7 (5) 
77'<7 



respectively, with the Pauli matrices a. t\ in Hj is de- 
fined by the pseudo-spin operator as 

rl = cos (^)T iz+ sin (^T ix , (6) 

with (m Xl m y ,m z ) = (1,2,3) where l(— x,y,z) indicates 
a direction of the bond connecting site i and site j. The 
magnitudes of the superexchange interactions J\ and J2 
satisfy the condition J\ > Ji > 0. The third and fourth 
terms in the Hamiltonian Eq. (JTJ describe the Hund cou- 
pling Jh{> 0) between the e g and t2 g spins 

T~Ch = —Jh / ] Si ■ Su, (7) 

i 

and the AF superexchange interaction Jaf{> 0) between 
the nearest neighboring ti g spins 

TLaf = Jaf / ] Su ■ Stj, (8) 

(ij) 

respectively. Su is the spin operator for the ti g electrons 
with S = 3/2. It was confirmed that this model is suc- 
cessful in description the spin and orbital ordered state 
in hole doped and undoped manganitesi^iiS 

Here, we examine the orbital ordered and disordered 
states in the unified fashion by applying the generalized 
slave boson method. The electron annihilation operator 
dija is decomposed into the holon hi, spinon s, ff and 
pseudo-spinon operators as 

with the local constraints of 

hthi + J2tl^ = l, (10) 

7 

and 

IX fl * = IX**r' (11) 

cr 7 

at each Mn site. We choose that both the holon 
and spinon are bosons, and the pseudo-spinon is a 
fermioni2&i2£i2£ This is because 1) in the ferromagnetic 
metallic phase, the hole carrier is recognized to be a 
fermion with orbital degree of freedom which is treated 
by h\ti^, and 2) the fermionic operator for the orbital de- 
gree of freedom ij 7 is suitable for describing the orbital 
disordered state. We adopt the mean field approximation 
where we take the bond order parameters for orbital 

E(*wi*W7')=*oxi, (12) 
77' 

and those for spin 
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and the site diagonal order parameters for orbital 

£(4r<'**r)=m?(0. (14) 
77' 

for /i = x and z, and those for spin 

="*;(*)■ (15) 

to in Eq. I|12[) is the electron transfer intensity between 
the nearest neighboring d 3z 2_ r 2 orbitals along the z axis. 
We also introduce a mean field for the holon 



(h\h i ) = -(hlh j )=x, 



(16) 



with the hole concentration x. In this scheme, the local 
constraints are released into the global ones: 



(17) 
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FIG. 7: The spin and orbital phase diagram in the hole con- 
centration x and temperature T plane. The parameter val- 
ues are chosen to be J2/J1 = 0.25, Jaf/Ji = 0.014, and 
to/Ji — 2. F, AF, P and L indicate the ferromagnetic phase, 
the antiferromagnetic phase, the paramagnetic phase, and the 
orbital liquid phase, respectively. 



and 

ia if 

By considering the strong Hund coupling Jh, the direc- 
tions of Si and Su are assumed to be the same. The 
effects of the lattice distortion and its coupling to the 
e g electrons are not considered explicitly in the present 
model. This is because, in the moderately hole doped 
region of our interest, the Jahn- Teller distortion is signif- 
icantly reduced as shown in the previous section. 

The phase diagram is obtained numerically by mini- 
mizing the free energy with respect to the order param- 
eters [Fig. 0. The parameter values are chosen to be 
J2/J1 =0.25, J A f/Ji=0.0U, and t /Ji = 2. For the 
Fourier transform of the site diagonal order parameters 
m^(k) and mf(fc), we assume that mf(7T7rO), mf(000), 
TOg(007r) and rn^(OOO) are finite in the ordered phases. 
Below the orbital ordering temperature Too , the site di- 
agonal orbital order parameters m^(k) appear. Tc and 
Tjv are the Curie temperature and the Neel temperature 
for the 4-type AF order, respectively, where to^(OOO) 
and mj(007r) appear, respectively. At x = 0, the stag- 
gered type orbital order mf (7T7t0) occurs far above T/v 
as observed in LaMnC>3 and PrMnC>3. By taking into 
account the cooperative Jahn- Teller effects which is not 
included in the present calculation, Too mav increase 
from the calculated value. With doping of holes, Too 
rapidly decreases and disappears around x = 0.15(= x c ). 
Above x c , m^(k) becomes zero, but the bond order pa- 
rameters xj(fc) still remain finite. That is, the orbital 
liquid state is realized above the critical point x c . As for 
the magnetic ordering, the ^4-type AF phase is changed 
into the FM phase through the canted AF one. The x 
dependence of Tq/Tn is almost flat below x = 0.1 and 



slightly increases with increasing x. The ground state 
electronic phase are classified into (I) the (canted) A- 
type AF phase with the staggered OO corresponding to 
LaMn0 3 and PrMn0 3 , (II) the FM phase with the stag- 
gered OO which is discussed later in more detail, and 
(III) the FM phase with the orbital liquid state corre- 
sponding to the FM-M state in heavily hole doped man- 
ganites, such as Lai-^Sr^MnOa with x ~ 0.3. In a sc- 
ries of Pri-^Ca^MnOa, it is interpreted that this phase 
does not appear because the charge ordered AF insulat- 
ing phase is stable down to the region of x ~ 0.3 instead 
of this FM-M phase. 

Now we focus on the FM phase with the staggered OO. 
We interpret that this corresponds to the FM phase in 
Pri-zCa-zMnOs with 0.15 < x < 0.3 of our present in- 
terest. The OO in this FM phase is of the staggered 
type where the two kinds of orbital are characterized 
by the orbital mixing angles (9/ — 6), as well as the 
OO in the A- type AF phase. This angle is defined as 
\0) = cos ^\d 3z 2_ r 2) + sin ||d 2 .2_j / 2). This is consistent 
with the azimuthal angle scan of the RXS experiments in 
P r o.75Cao.25Mn03; the azimuthal angle ip is the rotation 
angle of the sample around the x-ray scattering vector 
ki — kf in the RXS experiments. By measuring the ip de- 
pendence of the RXS intensity I(p), the information for 
the symmetry of the orbital ordered states is deduced^ 
The experimentally observed I(<p) in Pro.7sCao.25Mn03 
is proportional to a function sin 2 ip which is consistent 
with the theoretical prediction for I(tp) in the OO state 
of the (6/ — 9) type. 22 ' 34 The value of 9 in the present cal- 
culation is about 7r. That is, the orbital wave functions 
are given by {\d 3z 2_ r 2) ± \d x 2_ y 2))/^/2 where the elec- 
tronic clouds are more elongated along the z direction in 
comparison with the d 3x 2_ r 2 and d 3y 2_ r 2 orbitals. It is 
noted that, with doping of holes, the saturated orbital 
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FIG. 8: The spin and orbital phase diagram as functions of 
to/Ji. The parameter values are chosen to be J2/J1 = 0.25 
and Jaf/Ji = 0.014. 



PrMn0 3 and T c in PCMO with x = 0.3 is about 30% of 
the T N . 

In order to examine the spin dynamics in the orbital 
ordered FM phase, we calculate the spin stiffness by ap- 
plying the linear spin wave theory to the model Hamilto- 
nian in Eq. QJ. The effective exchange interaction J e $ 
corresponding to the interaction J; along the direction 
/(= x, y, z) in the Heisenberg model — Yltij) J{ Si ' Sj is 
obtained as 



reff 



2S 2 



(19) 



D\ is a stiffness constant arising from the transfer term 
of the Hamiltonian Tit given by 

Dj = xtoxl, 



(20) 



and Df is the one from Tij and TLaf given by 



Di 



moment \m^\ gradually decreases, and finally disappears 
at x c . On the other hand, the bond variable \x\\ grows 
up from zero at x = 0. That is, the two FM interactions, 
i.e. the double exchange interaction based on the kinetic 
term 7i t and the supcrexchange FM one from Ttj, co- 
exist. This is consistent with the present experimental 
results in Pro.7sCao.25Mn03 that the magnitude of the 
Jahn- Teller distortion, related to the magnitude of the 
diagonal order parameters , is significantly reduced in 
comparison with that in PrMnC>3. Unfortunately, in this 
calculation, this FM phase is not insulating, because the 
translational symmetry is assumed in the mean field ap- 
proximation. However, we numerically confirm that the 
band width of charge carriers, which is proportional to 
the bond order parameter Xt> is much reduced in the or- 
bital ordered FM phase in comparison with that in the 
FM phase orbital liquid phase. We believe that doped 
carriers in such a narrow band tend to be localized by 
other ingredients which are not taken into account in the 
present model. 

We also examine the phase diagram by changing the 
relative ratio of the kinetic energy and the exchange en- 
ergy to/Ji(2) [Fig- El- With decreasing to/<A(2), the crit- 
ical hole concentration x c , where Too disappears, in- 
creases and the region of the orbital ordered FM phase 
is extended to the region with higher x. On the other 
hand, Tc in the FM metallic phase is reduced and Tjv 
and Tc become less sensitive to x. These calculated 
results may explain a difference of the phase diagrams 
of Lai^Sr^MnOa (LSMO), Lai-^Ca^MnOa (LCMO) 
and Pri-^Ca^MnOs (PCMO). The parameter value of 
to/Ju2) in the calculation is expected to decrease in this 
order of the materials. The FM-I phase is realized in the 
regions of 0.1 < x < 0.15 for LSMO, 0.12 < x < 0.22 for 
LCMO and 0.15 < x < 0.3 for PCMO; the FM-I phase 
is extended and shifts to the region with higher x. Tc 
in the FM-M phase in LSMO is almost twice of Tjy in 
LaMnOa. On the other hand, difference between Tjv in 



(21) 



with 



J x = --(Ji-3J 2 ){l-x) 2 

+ \{Ji + J2)(\mf-\mf 



4xf ) - J 2 m$2) 



and 



~\{.h-z.h){i-xy 

\{Ji + J2) (\mf-Axf 



2J 2 m z t . (23) 



Here, we define the average magnitude of the Mn spin as 
S = (1 — x)/2 + 3/2. The x dependence of the effective 
exchange interactions, J x (~ J y ) and J 2 , are presented 
in Fig. 03 In the ^4-type AF phase near x = 0, J z is 
negative and its magnitude is much weaker than that of 
\J X \. This is experimentally confirmed by the inelastic 
neutron scattering in LaMnOa. With doping of holes, 
the sign of J z is changed into a positive one, and J x de- 
creases and J z increases. That is, the anisotropy of the 
exchange interactions is much reduced. This tendency 
of the x dependence of the anisotropy of is consis- 
tent with the present neutron scattering experiments in 
Pri-^CaajMnOs, although, in the calculation, the weak 
anisotropy of J e $ remains in the orbital ordered FM 
phase. Such an anisotropy may not be detected due to 
the coarse resolution of the present neutron scattering 
experiments. 



IV. CONCLUSION 

We have systematically investigated the nature of the 
FM-I state of perovskite manganites by a neutron scat- 
tering study on a single crystal of Pro.75Cao.2sMn03 and 
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FIG. 9: The x dependence of the effective exchange interac- 
tion and JtF . The parameter values are chosen to be 
J 3 /Ji = 0.25 and Jaf/Ji = 0.014. 



stiffness is similar to those observed in other FM-I man- 
ganites and has a smaller value than that in the FM- 
M phase, indicating the ubiquitous origin of the FM-I 
state. We theoretically construct the comprehensive spin 
and orbital phase diagram for the FM phase, which can 
explain the difference of the phase diagrams of several 
manganites. 

Lastly, we should note the importance of our work 
for the research of the I-M transition or the CMR phe- 
nomenon in hole-doped perovskitc manganites. Many 
previous experimental researches focusing attention on 
the phase separation as an origin of the I-M transition 
attribute an isotropic phase or a FM phase to a metal- 
lic phase. However, the existence of an isotropic ferro- 
magnetic insulating phase means one must pay better 
attention to identifying the metallic phase and consider 
the effect of the staggered type OO especially in the low- 
doped region. 



theoretical calculations. We have found that 1) the in- 
sulating behavior is robust against an external magnetic 
field, 2) the FM moment is almost fully polarized, 3) the 
Jahn- Teller lattice distortions are strongly suppressed, 
and 4) the spin exchange interactions are isotropic. These 
features are against interpretations of insulating behav- 
ior based on COs or phase separation, and are consistent 
to the staggered type orbital ordered state. The spin 
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